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ABSTRACT 

In this concluding part of the series of three papers dedicated to the 
Swift/B AT hard X-ray survey (BXS), we focus on the X-ray spectral analysis 
and statistical properties of the source sample. Using a dedicated method to ex- 
tract time-averaged spectra of BAT sources we show that Galactic sources have, 
generally, softer spectra than extragalactic objects and that Seyfert 2 galaxies are 
harder than Seyfert Is. The averaged spectrum of all Seyfert galaxies is consistent 
with a power-law with photon index of 2.00±0.07. The cumulative fiux-number 
relation for the extragalactic sources in the 14-170 keV band is best described by 
a power-law with a slope a = 1.55 ± 0.20 and a normalization of 9.6±1.9 x 10^'^ 
AGN deg"2 (or 396±80 AGN all-sky) above a flux level of 2xl0""erg cm^^ s"^ 
(~0.85 mCrab). The integration of the cumulative flux per unit area indicates 
that BAT resolves 1-2% of the X-ray background emission in the 14-170 keV band. 
A sub-sample of 24 extragalactic sources above the 4.5 a detection limit is used 
to study the statistical properties of AGN. This sample comprises local Seyfert 
galaxies (z=0.026, median value) and ~10% blazars. We find that 55% of the 
Seyfert galaxies are absorbed by column densities of Nh > lO^^H-atoms cm~^, 
but that none is a bona fide Compton-thick. This study shows the capabilities 
of BAT to probe the hard X-ray sky to the mCrab level. 
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1. Introduction 



There is a gen eral consensus tliat t iie cosmic X-ray background (CXB), discovered more 
than 40 years ago (iGiacconi et al.lll962l ). is produced by integrated emission of Active Galac- 
tic Nuclei (AGN). Popul ation synthesis m odels have successfully shown, in the context of 
the AGN unified theory jAntonuccil[l993h . that AGN with various level of obscuration and 



at different redshifts account for 80- 100% of the CXB below 4keV ( IComastri et al.l Il995 



Gilli et al.ll200ll : lTreister &: Urryll2005l ). Notwithstanding all the advances in the field a major 
question remains, do Compton-thick s ources exist in the numbers that seem to be required 
by population synthesis models (e.g. IComastri et al.l Il995l : iGilli et al.l l200ll ) to reproduce 
the shape of the CXB emission? Indication of the existence of s uch population comes from 
the analysis of the CXB fraction which is resolved-into-sources; IWorsley et al.l (120051 ) find 
that such fraction decreases with energy and that the unresolved component is consistent 
as being the emission of a yet undetected population of Compton-thick AGN. In synthesis, 
much evidence points towards the existence of Compton-thick AGN while only a handful of 
them are known and studied. 

The >10keV energy range is the most appropriate band for studying and selecting an 
unbiased (with respect to absorption) sample of AGN. This band is also the optimum band 
for the detection of Compton-thick objects. These elusive objects could have been missed 
because of the difficultie s of performing sensitiv e imaging of the hard X-ray s ky. The Burst 



Alert Telescope (BAT) (IBarthelmy et al.ll2005l ). on board the Swift mission ( IGehrels et al. 



20041 ) represents a major iraprovem ents in sensitivity for X-ray imaging of the hard X-ray 
sky. We refer to lAjello et al.l (120071 ) for details about the BXS survey. 



We applied an innovative image reconstruction algorithm to 8 months of survey BAT data; 
our survey covers ~7000deg^ reaching a limiting sensitivity of < 0.9 mCrab. This makes 
it one of the most sensitive survey ever performed in the hard X-ray domain. We detected 
49 hard X-ray sources of which 37 were previously unknown as hard X-ray emitters. Cor- 
relation with X-ray catalogs allowed us to identify 15 sources, while Swift/XBT pointed 
observations provided identification for another 15 obje cts. Furthermore, we optically iden- 
tified three new extragalactic sources (IRau et al.ll2007l ). Here, we investigate the spectral 
and statistical properties of all complete source sample. 



The paper is organized as follows. In section [2] we present the X-ray spectral analysis of 
the BAT sources. The details of the dedicated spectral extraction method are presented in 
the Appendix IA.4I We use the source spectra to build an X-ray color-color plot which is 
used to understand the mean properties of the source populations. In section [3l we apply 
the V/Vmax method to test the completeness of the extragalactic sample which is then used 
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to derive the number-flux relation. The section ends with a discussion about the statisti- 
cal properties of the extragalactic sample. Finally, we discuss the BAT results in section 
lU Throughout this work we use Hq = 70km s~^ Mpc~^ (hyo =1), k = 0, fimatter=0.3 and 
Ao=0.7 and the luminosities are given in erg s~^ h~^. 



2. Spectral analysis 

We have developed a dedicated spectral extraction method which allows to derive the 
time-averaged spectrum of all sources. The reader interested in the method is referred to 
the Appendix |X] for details. Using this method, we derived for all our source candidates 
a 6 channels energy spectrum in the 14-195 keV. The energy channels used in this analysis 
are (in keV): 14-22, 22-30, 30-47, 47-71, 71-121, 121-195. The energy bins were optimally 
chosen to produce similar error bars (in the different energy bins) for sources with power-law 
spectra. 

We found that 21 sources had at least soft X-ray observations by Swift/XKT or ASCA. For 
these sources, we jointly fit XRT/ASCA and BAT data. When fitting a source spectrum, we 
have preferred the simplest model yielding a good description of the data. The normalization 
of the ASCA spectra was allowed to vary (with respect to the BAT ones) to cope with the 
different epochs of the observations. This was not required when fitting XRT and BAT data. 
In general, the BAT spectrum of Galactic sources is well fitted by a thermal bremsstrahlung 
model. Instead, AGN are usually better described by a single power-law model. However, 
when <10keV data were available the fit required additional components (i.e. black body 
component for soft excess and/or gaussian model for the iron line). The detailed analysis is 
reported in Appendix IA.4I while the spectral parameters are summarized in Table [U 

The properties of the source sample can be studied using hardness ratios. We have, 
thus, defined HRi and HR2 as: 

HR ^^dium — hard soft — medium 

medium + hard ' soft + medium ' 

where the soft, medium and hard bands are respectively (in keV): 14-30, 30-71, 71- 
195. The hardness ratios, shown in FigH], are normalized to the range -1 and +1. Different 
symbols indicate different source classes. We also indicate the loci occupied by sources with 
a power law index in the range 1.0-3.5, or a bremsstrahlung spectrum with a temperature 
of 5-50 keV. A few things can be derived by the study of the hardness ratios. Galactic 
sources, which are usually characterized by soft X-ray spectra, have HR2 values < —0.3 
and HRi < —0.5 which is the typical region for sources with steep photon index. Indeed, 
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the five Cataclysmic Variables (CVs) present in the sample are all well fit by a relativistic 
bremsstrahlung model with a mean plasma temperature of 23 keV. 

Similarly we note that Sy2 galaxies seem to have (given the large uncertainties) harder X- 
ray spectra than Syls (larger values of HRi). The fact that type-2 AGN have systematically 
harder spectra than type-1 could be an evidence of the intrinsic difference between these 
two classes of objects. In order to study this in more detail we performed a stacked spectral 
analysi^ grouping the Seyfert galaxies detected by BAT into three classes: Seyfert 1, 
Seyfert 2 and intermediate Seyferts. The results, which are summarized in Table Ej show 
that the mean photon index of Seyfert- 1 and Seyfert-2 are different at mor e than 2 a level 



The same trend was previou sly noted in Seyfert gal axi es detected by OSSE f Zdziarski et al. 



20001) and by INTEGRAL flBeckmann et all 120061 ) . IZdziarski et al.l (120001 ) find that the 
difference in spectral index could be due to the different viewing angle between Seyfert Is 
and 2s. Indeed, the strength of Compton reflection decreases with the increasing viewing 
angle. Since the spectrum from Compton reflection peaks at 30keV followed by a steep 
decline, the larger the reflection component, the softer the spectrum. We tested this scenario 
using for Seyfert Is (Seyfert 2 s are successfully fitted by a simple power-law) the PEXRAV 
(IMagdziarz fc Zdziarskilll995l ) model in XSPEC. Indeed we get a good fit (x^ = 1-2/3) with 
a (minimum) reflection strength R> 1.1 (upper limi t is un cons trained by the fit) which is i n 
good agreement with what found by lZdziarski et al.l (120001 ) and lDeluit fc Courvoisierl (120031 ). 
Thus, the BAT data seem to confirm the larger reflection component present in Seyfert 1 
galaxies (with respect to Seyfert 2s) in agreement with the AGN unified model. Even though, 
the reflection component improves the fit, it does not affect the photon index of Seyfert Is 
which remains 2.30±0.12. 

We note, however, that most of the Seyfert Is (6 out of 9) have a low value of HRi 
denoting a steep spectrum. We thus tried to fit the stacked spectrum with a cutoff power- 
law model of the form Yj'^ e~'^^/^''\ Since the power-law index and the e-folding energy E^ are 
highly correlated we fixed the photon index to 2.0 (see below). The best fit e-folding energy 
is llO.S^'^glo keV (90% C.L.) with a reduced chi-squared which is substantially better than 
the one of the power-law model (0.9 vs. 1.4 and F-test probability of 0.08). The presence 
of a cutoff at ~100 keV in the X - ray sp ectra of Seyfert 1 seems also to be confirmed by the 
analysis of iDeluit fc Courvoisierl (120031 ) . 



Finally, we performed the stacked spectral analysis of all the Seyferts to investigate the 
averaged spectrum of the local AGN detected by BAT. The stacked spectrum, shown in 
Fig. [21 is consistent (in the 15-200 keV range) with a power-law model of photon index of 



^AU stacked spectral analysis are performed doing a weighted average of the spectra. 
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2.00±0.07 (90% C.L.). 



3. The hard X-ray extragalactic sample 



The extr agalactic sample, s hown in Table [3l was derived from the catalog reported 
in Table 2 of lAjello et al.l (120071 ) considering only objects at |6| > 15° and not spatially 
associated with the Large Magellanic Cloud. Here we describe the main properties of the 
sample. 



3.1. Completeness of the sample 



In order to compute the AGN number-fiux relation it is necessary to have a complete 
and un biased sample . Sinc e different regions of the sky have different exposure times, we ap- 
plied in lAjello et al.l (120071 ) a significance limit rather than a flux limit to define our sample. 
Now we want to test our extragalactic sample for completeness (i.e. derive the significance 
limit whi ch ensures to i nclude all objects above a given fiux limit) and we use the V/Vmax 
method ( Schmidt 1968 1^. This method which is applied to samples complete to a well- 
defined significance limit, can also be used to test the completeness level of a sample as a 
function of significance. For a significance limit below the true completeness level limit of 
the sample, the V/Vmax returns a value less than (V /'VMAx)true which would be the true 
test result for a complete sample. Above the completeness limit the (V/Vmax) values should 
be distributed around (V /V max) true within the statistical uncertainties. 



V/Vmax is computed for each source as {F/(atestSF))-^/^, where F is the fiux, SF is the 
la statistical uncertainty, atest is the significance level tested for completeness (and thus 
the term atest^F is the limiting fiux of the sky region where we detected the source), and 
the exponent —3/2 comes from the assumption of no evolution and uniform distribution 
in the local universe. (V/Vmax) is computed as an average of all sources detected with 
S/N> atest- For a gi ven mean value m = ( V/Vmax) and n sources, the error on (V/Vmax) 
can be computed as (lAvni fc Bahcalllll98C 



1/3 — m — 



n 



(2) 



^ In this test V stands for the volume where the object has been detected and Vmax is the accessible 
volume in which the object, due to the flux limit of the survey, could have been found. In case of no evolution 
<V/Vmax> = 0.5 is expected. 
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The results of the test are shown in Fig. [21 We find a constant value for significances 
> 4.5(7. The deviation from the expected 0.5 value is insignificant being less than lajfl 

We also remark that for completeness we are referring to the threshold above which all 
sources above the corresponding fiux limit are included in the sample. Furthermore, given 
the small redshift of the sample (see section the hypothesis of no evolution is justified. 



3.2. Extragalactic source counts 

The cumulative source number density can be computed as: 



Ns . 

i=i * 



(3) 



where A^s is the total number of detected sources in the field with fiuxes greater than 
S and flj is the sk y coverage associated to the flux of the i*'' source (shown in Fig. 9 of 
Ajello et all jiooj). 

The cumulative distribution is reported in Fig. HI We performed a maximum likelihood fit 
to the cumulative counts assuming a simple power-law model of the form A^(> S) = AS^"'. 
Here, A is the normalization at 2 x 10~^^ erg c ni~^ s~^ and a is the slope. As conventional 



we used the maximum likelihood estimator (e.g. I Crawford et al.lll97Cll ) to determine the best 
fit values. The normalization is not a parameter of the fit, but is obtained assuming that the 
number of expected sources from the best fit model is equal to the total observed number. 
The Poissonian error on the total number of sources provides a reliable estimate of its error. 

In the 14-170 keV band the best fit parameter is: a = 1.55 ± 0.20 (with normalization 
9.6±1.9 X 10~^ deg~^). The source count distribution is thus consistent with a pure Euclidean 
function {a = 3/2). From our data we expect that the number of all-sky AGN is 396±80 
brighter than 2 x 10^^^ erg cm~^ s^^. This corresponds to an integrated fiux of ~ 5 x 10^^^ 
erg cm~^ s~^ deg~^ or ~1.5% of the i ntensity of the X-r ay background in the 14-170 keV 



energy band as measured by HEAO-1 (iGruber et al.lll999l ). 



•^It is not uncomm on for coded mask detectors to produce a test value slightly above 0.5(e.g. 
Beckmann et al.l 120061 ). This is likely caused by systematic errors which tend to increase the (V/Vmax) 
value. 
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We can compare the surface density of extragalactic objects found by BAT with previous 
measurements by converting the BAT fluxes to other energy bands assuming a power law 
spectrum with photon index of 2.0 (see Section [2]) and evaluating the surface density above 
10~^^erg cm~^ s~^. 

The results of such comparisons are shown in Tab HI The BAT surface density is in 
agreement with the reported measurements, except for the case of the 0.5-2 and 2-10 keV 
surveys. Indeed, such surveys, at limiting fluxes of 10~^^ erg cm~^ s~\ are biased against the 
detection of absorbed sourceqj. It is also w orth noting that the re cent XMM measurement 
of the 5-10 keV source counts distribution (ICappelluti et al.l 120071 ) is in perfect agreement 
with our estimate. 



3.3. Statistical properties 

Above the 4.5 a the extragalactic sample, shown in Table[31 contains 24 AGN. 19 objects 
are classified as Seyfert galaxies, 3 as blazars, 1 as X-ray bright Optically Normal Galaxy 
(XBONG) and 1 as Quasar. The identification completeness of such sample is thus 100%. 

Excluding the blazars, the median redshift of the sample is z=0.026 (mean is z=0.046) 
giving a median luminosity of 10^^'^ erg s~^ (mean is 10^^'^ erg s~^) in the 14-170 keV band. 
Assuming a hydrogen column density of 10^^ atoms cm~^ as the threshold between absorbed 
and unabsorbed objects, we find that intrinsic absorption is present in ~55% of the sample. 
This fraction is lower than the 75% expected by the standard u nified model, which is derived 



by the opening angle of ionizations cones (e.g. lEvans et al.lll99ll ). However, this unexpectedly 



low fraction of absorbed AGN in the local Universe does not see m to pose any parti cular 



problem for the understanding and the synthesis of the CXB (e.g. ISazonov et al.l 120071 ) 



In Fig. [5]we show the intrinsic column density of the sources as a function of unabsorbed 
luminosity in the BAT band. Excluding the lower limits on the absorption, we do not find 
evidences of an anticorrelation between luminosity and absorption. We also note the presence 
of a rare very luminous (L^^ ~ 10^^ erg cm~^) highly absorbed (J^^h ~ 10^'^atoms cm~^) type- 
2 QSO. If the lower limits on the absorption will be confirmed, the total fraction of such 
objects might be in the range 5-15%. 

None of the source in Tab. [31 having 2-10 keV measurement, is a Compton-thick AGN. 
Our claim is supported by several evidences: 



''The bias decreases in deep fields and thus at lower fluxes (and higher redshifts) because the photoelectric 
cut-ofF is redshifted at lower energies. 
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As shown by lMatt et al.l (119971 ). for the case of NGC 1068, the spectra of Compton-thick 
AGN might be reflection-dominated (i.e. the reflection component is larger than the 
trasnmitted one). We thus tried to fit to each source a pure reflection model (PEXRAV 
in Xspec). For all the sources, except Mrk 704, th e fit is statis t ically unacceptable. 
However, Mrk 704 is not a Compton-thick source as iLandi et al.l (120071 ) have recently 
shown. 

Compton-thick source s generally show iron lines with equivalent widths of ~1 keV (e.g. 
Guainazzi et al.ll2005l ). The spectral analysis (see also values in Tab. [3]) shows that all 
sources have iron line equivalent widths smaller than 1 keV. 



The thickness parameter T, defined as Lg^j^Qj^^y/ Lqjjj (see also lBassani et al.lll999l ). 
can be used to identify Compton-thick sources (characterized by T<1 ). We coniputed 
the thickness parameter for all sources having OIII flux measurements ( iRau et al.ll2007l ) 
and 2-10 keV observations (see Tab. [3]). All sources except NGC 2992 (which is however 
unabsorbed) have thickness parameter values consistent with the values expected for 
Compton-thin AGN. 



We evaluated the radio-loudness of AGN using the R-index defined in iLaorl (120001 ) as 
R= f,y[5 GHz)//j,(4400A); the distribution of R-values has been shown to be bimodal with 
a minimum at R=10, commonly used to define radio-loud (above 10) versus radio-quiet ob- 
jects. Interestingly, we note that a relevant fraction (~ 40%) of the BAT AGN is radio-loud 
and that these objects show a systematically harder X-ray spectra than Seyfert galaxies 
(mean of 1.66 vs. 2.00). There is large co nsensus that radio-loud quasars host more mas- 
sive b lack holes than radio-quiet ones (e.g. iMetcalf fc Magliocchettill2006l : iMcLure fc Jarvis 
20041). Howeyer, th ere is no simple explanation for this radio-loudness dichotomy. Recently 



Sikora et al.l (120071 ) showed that the radio-loudness parameter inversely correlates with the 



Eddington ratio (fraction of bolometric to Eddington luminosity) for both spiral/disk and 
elliptical galaxies. The fact that spiral-hosted AGN are radio-quiet at high accretion lumi- 
nosities supports th e idea that the black hole spin plays a major role in the jet production 
dSikora et all [2OO7I1R. As a confirmation, we find a good correlation of intrinsic X-ray 
luminosity and radio-loudness (Spearman rank test being 0.57 with probability of 0.003). 
Such cor relation is expected if there is a fundamental connection between accretion and jet 
activity jMerloni et al.lboosh . 



^ In fact, by merging processes, black holes in elliptical galaxies are expected to have larger spins than 
those in spiral/disk galaxies. 



- 9 - 



Discussion 



We have used the BAT X-ray survey to study key properties of the local (z<0.1) AGN 
population. Our survey is based on the 14-170 keV fluxes and it is sensitive to AGN with 
column densities up to Nh ~ 5 x 10^^ atoms cm~^. Indeed, for a typical source with photon 
index of 2, the decrease in flux for column densities of Nh ~ 10^^ atoms cm~^ is only ~ 7% 
and ~ 55^ for column densities of ISIh ~ 3 x 10^^ atoms cm~^ . Thus, we can affirm that 
this survey is relatively unbiased with respect to photoelectric absorption. 



Most of the population synthesis models (lUeda et al.ll2003l : lTreister fc Urryll2005l : iGilli et al 



20071 ) predict that Compton-thick AG N (logN^ >24) provid e a significant contribution to 



the bulk of the CXB emission at 30 keV jMarshall et aL 



19801) ■ Although studies of the local 



Universe (e.g. Risaliti et al., 1999) have shown that Compton-thick objects should be as nu- 
merous as moderately obscured AGNs (logN// <24) and thus roughly 1 /3 of the tota l AGN 
population, only a handful of these sources are known (jComastril 120041 ) . iGilli et al.l (120071 ) 
estimate that the expected fraction of Compton-thick objects at limiting fiuxes probed by 



BAT and INTEGRAL (~ 10 erg cm s ) is in the 15-20% range. However the measured 
fraction of detected Compton-thick objects by these instr uments is, so far, close to, or less 
than, 10% jMarkwardt et ahlboosl : iBeckmann et allbooeh . 



The BAT extragalactic sample contains only one source SWIFT J0823. 4-0457, which 
given its colors (see Fig. [1]) might be Compton-thick. However, the joint XRT and BAT 
spectra show that the absorption is below the Compton-thick level (Nh ~ 10^^ atoms cm^^). 
We must therefore conclude that no Compton-thick AGN are present in our extragalactic 
sample. The probability of not detecting Compton-thick objects in a sample of 24 AGN 
when the expected fraction is 20% (15%) is ~0.007 (~0.03) while it is 0.1 if the expected 
fraction is 10%. These probabilities increase (0.03, 0.09 and 0.2 for the 20%, 15% and 10% 
cases) if we assume that the only source which lacks < lOkeV measurement (J0854. 7-1-1502) 
is Compton-thick. Thus, the BAT data discard tat >2 a level the hypothesis that Compton- 
thick AGN may represent a fraction of ~ 20% of the total AGN population. 

We find that Sy2s have harder spectra than Syl s in agreement with what has been de- 
duced from OSSE, BeppoSAX and INTEGRAL data ( Zdziarski et al .112000 : Deluit &: Courvoisier 
20031 : iBeckmann et al.ll2006l . respectively). We tested whether this difference could be ac- 
counted for by Compton reflection and/or by a high-energy cut-off. We flnd that the reflec- 
tion component improves the flt to the Syl averaged spectrum (the F-test shows that the 
reflection is signiflcant at more than the 92% level), but it leaves unaltered the photon index. 



'Photoelectric absorption as well as Compton scattering has been taken into account in this estimate. 
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Thus, the difference in photon indices among Syls and Sy2s cannot be ascribed solely to 
orientation effects (a stronger reffection is expected for face-on objects) . The spectra of Syl s 
show hints of a spectral cut-off at ~100keV in agreement with lDeluit &: Courvoisierl (120031 ). 
According to thermal Compton models, the absence of a cut-off in Sy2s might indicate a 
higher temperature of the Comptonizing medium (with respect to Syls) or that non-thermal 
Compton scattering plays an important role. Nevertheless, given the low S/N of our sources 
evidences for the cut-off in the Syl spectra are weak. 

The best power-law fit to the extragalactic source counts distribution yields a slope 
of a = 1.55 ± 0.20 which is consistent with an Euclidean distribution. From the best 
fit, we derive a surface density of AGN of 9.6 ± 1.9 x 10~^deg~^ above the flux limit of 
2xl0~^^erg cm~^ s~^; this estimate is in very good agreement, when converted to the 20- 
40ke V band, with the rece n tly derived source coun ts distribution based on INTEGRAL 
data (IBeckmann et al.ll2006l ). iBeckmann et al.l (120061 ) find a slope of 1.66±0.11 which is also 
consistent with our measurement, but steeper than the 1.5 Euclidean value. Even though 
this could be due to a non-perfectly computed sky coverage, the authors suggest that the 
distribution of AGNs in the local Universe may not be isotropic because of local clustering 
of sources (e.g. the local group of galaxies). 

The BAT source count distribution resolves only 1-2% of the CXB into extragalactic 
sources; nevertheless as it is unbiased with respect to absorption it gives important informa- 
tion relative to the fraction of obscured sources which are missed by deep < lOkeV surveys 
because of absorption. The extrapolation of the BAT source count distribution to the 2- 
10 keV band assuming an unabsorbed spectrum with photon index 2 yields a surface density 
of AGN of 1.6±0.32 x 10~^deg~^ above 10~^^erg cm~^ s~^; wh ile the surface density as 



extrapola ted to brighter fl uxes by XMM (ICappelluti et al.l 120071 ) and as predicted by the 
model of iGilli et al.l (120071 ) is 0.9x lO^^deg"^. The factor ~ 2 more sources BAT sees can 
be explained in term of absorptio n. Indeed, if we t a ke int o account the absorption distri- 
bution derived for BAT AGNs by iMarkwardt et al. J2005h (thus assuming that 66% of all 
AGN are absorbed with a mean column density of 10^^ atoms cm~^) we get a surface den- 
sity 0.86±0.17 X 10~^deg~^ which is consistent with the XMM extrapolation and the model 
prediction. 

The extragalactic sample is composed of ~90% emission-line galaxies and ~10% blazars. 
We find that 55% of the emission-line galaxies are obscured by absorbing columns larger than 
10^2 H-atoms cm 2. This fraction is in agreement with the INTEGRAL measurements (e.g 



Sazonov et al. 



20071). bu t less t han what is suggested (~ 75%) by the unified AGN model. 



However, ISazonov et al.l (l2007l ) successfully showed that low-luminosity (mostly absorbed) 
AGN account for much, ~90%, of the luminosity density of the local Universe. This is 



- 11 - 



also confirmed by the model of iGilli et al.l (120071 ) which shows that the required fraction of 
obscured sources varies with intrinsic luminosity being 3.7 and 1.0 below and above lO^'^'^ erg 
s~^. A relevant fraction (~40%) of the BAT-detected AGN is radio- loud. These objects 
show a systematically harder X-ray spectra than Seyfert galaxies (1.66 vs. 2.00). The hard 
photon index and the correlation of radio-loudness with X-ray luminosity suggest that a jet 
is presently at work in all these objects. Our sample also comprises 1 (and possibly up to 
3 considering the ROSAT lower limits on the absorption) highly luminous highly absorbed 
QSO. 



5. Summary 

We use the Swift/BAT instrument to study the properties of the local (z<l) AGN in 
connection with the synthesis of the X-ray background emission. The results of this study 
can be summarized as follows: 

• Despite the consensus that Com pton-thick objects rnay represent a su bstantial fraction 



of the local AGN population (e.g. lRisaliti et al.lll999l : lGilli et al.ll2007l ). we do not detect 



any such object. The probability associated to this non-detection is 0.007, 0.03 and 
0.1 when assuming that their fraction should be 20%, 15% and 10% of the total AGNs. 
BAT discards at > 2 cr the hypothesis that the fraction of Compton-thick objects is 
20%. 

Seyfert 2 galaxies have harder X-ray spectra than Seyfert 1. We find that this difference 
cannot be ascribed solely to the different viewing angle and thus to the different amount 
of Compton reflection which is expected. The Seyfert 1 galaxies comprised in our 
sample show weak evidences for a spectral cut-off in the ~100keV range. This might 
highlight an intrinsic difference among the two classes. Indeed, the absence of a cut- 
off in the spectra of Seyfert 2s might indicate a different (higher) temperature of the 
Comptonizing medium or that non-thermal Compton scattering play an important 
role. 

The best power-law fit to the extragalactic source counts is consistent with a Euclidean 
function with slope of 1.55±0.20. At the current limiting fluxes (2xl0~^^ erg cm~^ s~^), 
BAT resolves only 1-2% of the CXB emission in the 14-170 keV band. 

The fraction of emission-line AGN which is absorbed by Nj;^ > 10^^ atoms cm~^ is 
~55%. This is lower than the 75% expected by the standard AGN unified model. 
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This work shows the capabihties of BAT to produce an unbiased sample of AGN which 
is important for the understanding of the synthesis of the CXB emission in the hard X-ray 
band. 
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A. Spectral extraction method 

We have developed a method to extract the averaged long term spectrum of a source. 
In this method, the spectrum is obtained as a weighted average of the source spectra of all 
observations where the source is in the field of view. In particular, the averaged source count 
rates in the i-th energy channel, Ri, and their error (Tj, are given by the following equations: 



(Al) 



where rj is the source count rate in the j-th observation, Wj is the weight used and the sums 
extend over all observations which contain the source. Using the inverse of the count rate 
variance Vj as a weight, the previous equations simplify to: 



However, the spectra entering in Eq. I A2 1 must be corrected for off-axis count rate varia- 
tion and for residual background contamination. We explain below the way these corrections 
are implemented. 



A.l. Rate variation as a function of off-axis angle 

The detected count rates strongly vary with the position of the source in the FOV; a 
source at the far edge of the partially coded FOV (PCFOV) can experience a decrease in 
rate of a factor 2 (depending also on energy) when compared to its on-axis rate. 

The standard Swift-BAT imaging software corrects for geometrical off-axis effects like 
cosine and partial coding (vignetting) effects; it is only when the response matrix is generated 
(tool batdrmgen) that other effects like detector thickness and effective area variation are 
taken into account. 

Since in equation IA2I we are averaging over spectra at different positions in the FOV, we 
need to take into account the variations in the rates produced by the detector response. In 
order to do so, we have analyzed a series of more than 1000 Crab observations. For each of 
our 6 energy channels we made a polynomial fit to the Crab rate as a function of the off-axis 
angle, and derived a set of corrective coefficients. These coefficients are then used to correct 
the rates of each source spectrum in order to transform them to the equivalent on-axis rates. 
The variation of the Crab rates as a function of position in the FOV is reported in Fig. [61 
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A. 2. Residual background contamination 



In order to extract a source spectrum from survey data (in form of Detector Plane His- 
tograms, DPH) the user must first produce a mask of weights (tool batmaskwtimg) for the 
source position and then use this mask to extract the detected counts from the array (tool 
batbinevt). The weights are chosen such that the resulting spectrum is already background- 
subtracted. This is an implementation of the standard mask weighting technique called 
balanced correlation ( iFenimore fc Cannon! Il978l ). The automatic background subtraction 
works as long as the noise in the array is fiat and not correlated with the mask pattern. 
These conditions are not always satisfied and a small background contamination can arise. 



The total background contamination for the case of the Crab is < 2% when compared 
to the Crab on-axis rate in the 14-195 keV band. Thus, this contamination does not pose 
problems for bright sources. However, it becomes relevant for the spectral analysis of faint 
objects with intensities of ~mCrab. 

In order to correct for this residual background contamination, we fit the batclean background 
model to each energy channel in order to create a background prediction for each of them. 
Convolving these background predictions with the mask of weights generated for the source 
under analysis yields the residual background term which the mask weighting technique did 
not manage to suppress. 



A. 3. Spectral fitting 

The final source rates in the i-th energy channel are computed as: 



rU''^ -b,)-K{E ,e)-l/V, 

Ef=o m 



= ^^=''' ' ' ' (A3) 



where bj is the residual background term, K{E, 9) is the parametrized instrumental 
response as function of the energy channel and the off-axis angle and Vj is the rate variance. 
The weighted aver aged spectrum is then input, together with a the BAT response matrix. 



to XSPEC 11.3.2 flArnaudlll996f ) for spectral fitting 



Finally, we checked that the averaged Crab spectrum, obtained with the above method, 
is consistent with the standard (BAT) Crab spectrum as detected in each observation (pho- 
ton index of 2.15 and normalization of 10.15 photons cm~^ s~^ at 1 keV in the 15-200 keV 
energy range). 
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A. 4. Notes on individual sources 



We report a brief description of tlie source spectra for all new, or interesting, sources 
found in tliis analysis. All quoted errors are 90%. The spectra of all the sources are reported 
in Figures [71[8l[9l[l0l[IIl[Il,[l3]and[Il (available online) . 

3C 105.0 is a Sy2 galaxy. The BAT and XRT data can be fitted by an absorbed power 
law model with photon index of 1.65±0.13 and hydrogen column density of 29.4^48 x 10^^ 
atoms cm~^. Given its absorption and its luminosity (4.45x10^^ erg s^^), 3C 105.0 is a 
highly-absorbed highly-luminous QSO. 



1 AXG J042556-5711 (also known as IH 0419-577LB 1727, lES 0425-573 and IRAS 
F04250-5718) is a radio-quiet Seyfert galaxy which has been observed over re cent years by 
ASCA, ROSAT, BeppoSAX and recently also by XTE jRevnivtsev et allbooeh . The ASCA 
and BAT data are well fit by an unabsorbed cut-off power-law model with photon index of 
1.54±0.028 and cut-off at rSttHkeV. 



3C 120 is a Syl galaxy. This source was observed by ASCA. The best fit to ASCA and 
BAT data is an absorbed power law model with absorption consistent with the galactic one 
and photon index 1.80to;Q4 and a black body component with a temperature of 0.27to;o25 keV. 



MCG -01-13-0 25 is a Svl.2 (in NED, but Syl in SIMBAD) galaxy detected in soft X- 
rays by ROSAT fjVoges et al.lll999l ). The BAT spectrum is consistent with a power law with 
photon index of 1.6l;o;47 and it extends up to 200 keV. 



SWIFT J0505. 7-2348, also know as XSS J05054-2348 flRevnivtsev et al.ll2006f ) is a Sy2 

galaxy. When combining both XRT and BAT data for this source we get an intrinsic, rest 
frame, absorption of 4.8[!!q7 ^ ^O^^ atoms cm~^ and a photon index of 1.77^ 



-0.08 
-0.07- 



CSV 6150, also known as IRAS 05078+1626, is cataloged as Syl in SIMBAD and as Syl.5 
in NED. The BAT spectrum can be fitted with a power law with photon index of 1.941q'23. 
The source flux in the 14-170 keV band is 6.3~^^q x 10^^^ erg cm~^ s"^ while the luminosity 
is AAtli X 10^3 grg s-\ 



4U 513-40 is a low mass X-ray binary detected in X-rays by EXOSAT (iGiommi et al. 
19911 ). The BAT spectrum can be fit by a bremsstrahlung model with temperature of 
29.7^^1 keV. 
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QSO B0513-002 is a Syl galaxy. The BAT and ASCA spectra can be fitted by an absorbed 
power-law model and a black body component. The required absorption is in agreement with 
the galactic one. The photon index and the plasma temperature are respectively 1.83lo;o^g 



and 0. 27^002 keV. We also detect a an iron line whose equivalent width is OO.Si^ygy eV. 



+66.2 



SWIFT J0517.1+1633 is a new hard X-ray source (lAiello et all 120071 ). The BAT spec 



trum is best fit by a power law model with photon index of 2.0 



+0.23 
-0.26- 



ESO 362- G 018 is a Syl galaxy detected at hard X-ray by BAT flTueller et al.ll2005[ ). 
The BAT and XRT data are best fit by an absorbed power-law model with photon index of 
1-501;o;q2 and absorption consistent with the galactic value. 



Pictor A is a radio-loud Syl galaxy initially detected in X-ray by the EINSTEIN ob- 
servatory (lElvis et al.lll992l ). The best fit to ASCA and BAT data is an absorbed power-law 
model with photon index of I.8I002 and intrinsic absorption of 1.141q'q5 x lO^^atoms cm~^ 
slightly in excess of the galactic one ( 4 x 10^° atoms cm~^). 



ESO 362-G021 is a BL Lac object. ASCA and XRT data are available for this source. 
The best fit to ASCA, BAT and XRT data is an absorbed power law with photon index of 



1.721;q q4 and intrinsic column density of 0.14;^q q2 x 10^^ atoms cm 



0.02 



TV Col is a DQ Her type cataclysmic variable already detected at soft and hard X- 
rays. A power law fit to the BAT spectrum does not yield acceptable results; instead a 
bremsstrahlung model with a plasma temperature of 28.2^3;8keV fits the data well. 



TW PIC is a cataclysmic variable of the DQ Her type (INorton et al.ll2000l ). The BAT spec- 



trum is best fitted by a bremsstrahlung model with a plasma temperature of IS.Slg^g^ keV. 
The fiux of 5.5xl0~^^erg cm~^ s~^ in the 20-40 keV band is a factor 2 lower than the one 



reported in a recent INTEGRAL measurement (IGotz et al.l 120061 ) suggesting variability. 



LMC X-3 is high mass X-ray binary (HXB). The BAT spectrum is consistent with a power 
law whose photon index is 2.0^ 



-0.4 
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LMC X-1 is a well known black hole candidate. It is detected up to 200 keV with a steep 



photon index of 2 .3ln 9n- The fiux is a factor 2 lower than the one measured by INTEGRAL 



( iGotz et al.ll2006l ). suggesting variability. 



PSR B0540-69.3 is a young rotation-powered pulsar recently detected up to 60 keV also by 



INTEGRAL f lGotz et alJbood : ISlowikowska et alJbood ). The Pulsar is detected in BAT up 
to 200 keV and it is spectrum can be modeled as a power law with photon index of l.SS^glg. 



PKS 0537-286 at z= 3.1 is one of the most luminous high-redshift quasar. Recognized 
first as a ra dio source (jBolton fc Butlerl 119751 ) it was discovered in X-rays by the Einstein 
observatory (jZamorani et al.lll98ll ) and then studied by ROSAT, ASCA and lately by XMM. 
The BAT detection in hard X-rays is the first to date, however there i s a claim that PKS 0537- 



286 b e the MeV counterpart of the EGRET source 3EG J053 1-2940 (ISowards-Emmerd et al. 



20041 ) ■ A joint spectral fit to XRT and BAT data reveals an exceptionally hard spectral slope 



r . or +0.06 
Ol i.OO_oo8- 



PKS 0548-322 is a well known blazar already detected in hard X-rays (see for exam- 
ple Donato et al., 2005). A joint spectrum of XRT and BAT data with an absorbed power 
law model yields a photon index of 1.8]'^g;o3 and an intrinsic absorption of 2.57!l^o;5 ^ 10^° 
atoms cm~^. 



NGC 2110 is a well known Sy2 galaxy. The BAT, ASCA and XRT data can be fit by 
an absorbed power law model (photon index of 1.621q'q} and intrinsic hydrogen column den- 
sity of 4.0lo;o7 ^ 10^^ atoms cm~^) with a soft excess which could be described as black 
body component with temperature of 0.471q;q2- We also detected an unresolved Fe Kq, of 
equivalent width of 118^^53 eV. 

LEDA 75476, also known as 3A 0557-383, EXO 055620-3820.2 and GTS B31.01, is a 
Syl galaxy. The BAT spectrum is consistent with a power law model with photon index of 
2.0 ± 0.4. The ASGA and BAT are well fit by an absorbed power law model with photon 
index of 1.741q q3 and intrinsic absorbing column density of 2.2^0^3 x lO^^atoms cm~^. A 
clear excess below 2 keV is detected in the ASGA data and this can be modeled as a black 
body component with a temperature of 0.281o;q5 keV. A Fe Kq line is also required by the 
fit (F-test yielding a probabihty of the line being spurious of 10~^) and its equivalent width 
is 0.132 keV. The reduced of the overall fit is 1.1. 

ESO 490-G 26 is a Syl. 2 galaxy. The joint XRT-BAT spectrum can be described as a 
power law with photon index of 1.901o'q4 and an intrinsic, in addition to the galactic, ab- 
sorption of 2.7!!^Q 05 ^ ^^^^ atoms cm~^. The flux and the luminosity in the 14-170 keV band 
are 3.6+};3 x lO^^^erg cm~^ s"^ and 4.7t3;g ^ 10^^ erg s"^ 



SWIFT J0727. 5-2406 has a spectrum consistent with a power law model with photon 
index of 1.53 ± 0.54. As already noted by iRau et al.l (120071 ) this BXS source is likely asso- 
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ciated with the nearby ROSAT source IRXS J072720.8-240629 and with the radio object 
NVSS J072721-240632. 



V441 Pup is a high mass X-ray binary where the companion was optically identified as 
a Be star. The BAT spectrum is very steep and it can either be fitted by a power law 
with a photon index of 4.5±1.5 or by a bremsstrahlung model with a plasma temperature 
of 12.41^3.^ keV. 

BG CMi is a well known intermediate polar. The BAT spectrum is consistent with a 
bremsstrahlung model with a plasma temperature of 31.3l^^4 2 keV. 



SWI FT J0732. 5-1331 was detected for the first time by BAT in hard X-rays (lAjello et al. 
20061 ). It was then identified as a new intermediate polar (Wheatley et al. 2006 and refer- 
ences therein). The BAT spectrum is consistent with a bremsstrahlung model with a plasma 
temperature of 33.2^^4;2 keV. 



SWIFT J0739. 6-3144 is a newl y discovered har d X-ray source flAiello et al.l 120071 ) . re- 
cently identified as a Sy2 galaxy (iRau et al.l 120071 ) . A simple power law fit to the BAT 
spectrum yields a photon index of 1.771q;43 • ^^^o estimated the lower limit on the ab- 
sorbing column density considering the undetection by ROSAT; this limit is ~ 2 x 10^^ atoms 
cm~ . The fiux and the luminosity in the 14-170 keV band are 2.?,t\l x 10"^^ erg cm s 
and ^.2t\l x lO^^g^g g-^ 



SWIFT J0743.0-2543 is a newly discovered hard X-ray source flAiello et all 120071 ). The 
BAT sp ectrum is consist ent with a power-law model with photon index of 1.78lQgg. As 
noted in lRau et al.l (120071 ) this BXS source is likely to be associated with the ROSAT source 
IRXS J074315.6-254545 and the galaxy LEDA 86073. 



IGR J07597-3842 is a source first detecte d by INTEGRAL in th e VELA region fiden Hartog et al 
20041 ). It was identified as being a Syl.2 (iMasetti et al.ll2006bl ). This source was also ob- 
served by XRT and when jointly fitting XRT and BAT data we get that the best fit is an 
absorbed power law with photon index of 1.81q;q7 ^"^^ column density of 5.81q;5 x lO^^atoms 
cm~^ consistent with the galactic foreground absorption. The source is thus unabsorbed. 



The fiux and the luminosity in the 14-170 keV band are 4.2 
X 10^^ erg s'^ 



+°:^ X 10 



—11 —2 —1 

erg cm s 



and 



UGC 4203 is a Sy2 galaxy. As already noted in lMatt et al.l (120031 ). this source shows tran- 
sitions between a refiection-dominated and a transmission-dominated spectrum. The ASCA 
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and B AT data can be successfully fit by a reflection model (PEXRAV. lMagdziarz &: Zdziarski 



19951 ) with photon index of 1.68±0.1 and a reflection normalization of 65.2_23 gy and a promi- 



nent iron line with equivalent width of O.Tto ekeV. A soft excess at energies < IkeV can be 
modeled as a black body component with a temperature of 0.3lQQ5keV. 
XRT data are also available for this source. However, the XRT spectrum has a lower quality 
than the ASCA one. In the XRT observation, the source is found in a transmission-dominated 
state; the best fit model is an absorbed reflection model (the reflection component is required 
by the BAT spectrum) with hydrogen column density of Nj|/=12.5^37 x 10^^ atoms cm~^, 
photon index of 2.0]l^o25 and reflection normalization of 2.12^^5. 



SWIFT J0811. 5+0937 is a new BXS source detected in lAiello et al.l (12007 ). The BAT 



spectrum is consistent with a power law with photon index of 2.2^'^Qg. iRau et al.l (120071 ) iden- 
tified RX J081132. 4+093403 as possible counterpart. Optical spectroscopy revealed that this 
source is a candidate X-ray Bright Optically Normal Galaxy (XBONG). If we extrapolate 
the BAT power law to the ROSAT-PSPC energy band (0.1-2.4 keV) we get no indication of 
intrinsic absorption. 



SWIFT J0823. 4-0457 is a source detected for the first time in hard X-rays by BAT ari d 
associated, during an XRT follow-up, with the galax y FAIRALL 0272 (lAiello et al.l 120071). 
An optical follow-up showed that the source is a Sy2 (IMasetti et al.ll2006al ). XRT and BAT 



data are best fit by a highly absorbed power law. The photon index is 1.84]'^q22 and the 
absorbing column density is 19.3]t5i x 10^^ atoms cm~^. 



Vela PSR has a spectrum consistent with a power law whose photon index is 1.88+0.2. 



FRL 1146 is a Syl galaxy deteced in hard X-ray by INTEGRAL (IBird et all 120061 ). The 
BAT spectrum is characterized by a power law with photon index of 1.88^*^03]' extending 
up to 200 keV. The 14-170 keV flux and luminosity of 3.3lo;7 x 10"" erg cm"^ s~^ and 
7.2lj^4 X 10^'^ erg s~^ are in agreement with the INTEGRAL measurement. FRL 1146 was 
also detected in the ROSAT all-sky survey at 12 count s~^, considering the extrapolation 
of the BAT power law to the ROSAT band yields ~8 count s~^ so it is very likely that the 
source in unabsorbed. 



3C 206 is a narrow line, radio loud, QSO detected for the first time in hard X-rays (>20 keV). 
It was detected by Lawson & Turner (1997) using GINGA in the 2-10 keV The BAT spec- 
trum is consistent with a pure power-law model with photon index of 1.9 5ln'qQ- 3C 206 was 



detected by the ROSAT PSPC with 0.37 ct s"^ during the all-sky survey (IVoees et al.lll999r ): 
if we use the BAT power law spectrum and we extrapolate it to the 0.1-2.4 keV band, we 
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get that no additional absorption (with respect to the galactic one) is required to match the 
observed ROSAT count rate. 



SWIFT J0844.9-3531 is a new hard X-ray source detected in lAjello et all (120071 ). The 

BAT spectrum is consistent with a power law model wit h photon index 1. 911° ^^. The flux 



in the 14-170 keV band is 1.7+o;g x 10"^^ erg cm'^ s" 
BXS source might likely be associated with the ROSAT source IRXS J084521. 7-353048. 



Rau et al.l (120071 ) noted that this 



SWIFT J 0854.7+1502 is a new hard X-ray source detected in lAiello et all (120071 ) and 
identified in Rau et al. (j2007 ) as a Sy2 galaxy. It has a flat spectrum which can be modeled 
as a power law with photon index 1.41^'^Q'g. A lower limit on the absorbing column density of 
5x 10^^ atoms cm~^ can be derived by the non-detection of this source in the ROSAT all-sky 
survey. 



SWIFT J0917.2-6221 is a new hard X-ray source. We analyzed a 7 ks XRT observa- 
tion of this source. The XRT and BAT data are well flt by an absorbed power law model 
whose photon index is 1.87lo;o4 absorbing column density of 1.33lo;|o ^ ^0^^ atoms cm~^. 
A clear excess is present at energies <lkeV and this can be well described as a black body 
component peaking at 0.14 keV. The flux and the luminosity in the 14-170 keV band are 

erg s~^. 



2.6+o:8 X 10"^^ erg cm'^ s 



1 and 20.0+^:° x 10^^, 



Mrk 0704, or S WIFT 0918.5-1-1618, is another source found thanks to our algorithm 
(lAjello et al.l 120071 ). During an XRT follow-up, the galaxy Mrk 704 was fou nd as the BAT 



counterpart. Mrk 704 was previously detected in soft X-rays by R OSAT (ISchwope et al. 



20001 ). In a recent optical follow-up, the galaxy was found to be a Syl (iMasetti et al.ll2006al ). 
We have analyzed ASCA, XRT and BAT data for this source. The best flt to the three 
datasets is a a partial covering model where the covering fraction is 0.5 and the powerlaw 
photon index is I.SG'^q I^. The source is highly absorbed with a column density of 1.5^q;3 x 10^^ 



atoms cm . We also detected an iron line whose equivalent width is 160 eV. 

4U 09 19-54, detected at very high signiflcance, is a LMXB also known to produce X-ray 
bursts (jjonker et al.ll200ll ). Its spectrum is characterized by a steep photon index of 2.35 ± 



0.25, alternatively a bremsstrahlung model with a plasma temperature of 45.11^9^80'^ keV 



yields a better x^- 



MCG -01-24-012 is a Sy2 galaxy already detected in hard X-rays by Beppo-SAX (IMalizia et al 
2OO2I ). When fltting both XRT and BAT data we get that the spectrum is consistent 



with an absorbed power law whose photon index is 1.7lo'o7 intrinsic absorption is 
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6.5lo:? X 10^2 atoms cm-^. 



NGC 2992 is a Sy 1.9. 



The best fit for combined XRT, ASCA and BAT data is an 

06 



absorbed power law witli plioton index of 1.24]']qq5 and intrinsic liydrogen column density 
of 0.17: 

-ao3 ^ 10^2 atoms cm . We also detected the presence of an unresolved Fe Kq, line 
whose equivalent width is 0.52lQ 'j' keV in agreement with an old Beppo-Sax measurement 
(IGilli et al.ll200ll ) where the reported column density is IxlO^^atoms cm"^. 



ESQ 434- G 040 is a known Sy2 galaxy recently detected in hard X-rays also by INTEGRAL 
( iBird et al.ll2006l ). A joint fit to ASCA, XRT and BAT data with an absorbed power law 
model yields a photon index of 1.77;^q;q7^ and a column density of 1.5lQ;Qg^ x 10^^ atoms cm"^. 
A clear excess below 2 keV can be modeled as a black body component with a temperature 
of 0.13;°;°J^. An iron K„ 
the line being spurious is ~ 10"^'' 



line, with an EQW=85.5t33, ^^^^ detected. The probability of 



3C 227 is a Syl galaxy and also a Radio galaxy. The BAT spectrum is consistent with 
a power law model of photon index l.QGtasg- This source wa s detected at a level of 0.0 16 ct 
s-^ in a 11 ks long ROSAT-PSPC observation (0.1-2.4 keV) jCrawford fc Fabian! [l995l ) . In 
order to match the ROSAT observed count rates, the extrapolation of the BAT power law to 
the 0.1-2.4keV band requires an absorbing column density of at least 5 x 10^^ atoms cm"^. 
A recent Chandra observation confirms that 3C 227 is indeed an absorbed Syl. However 
the joint Chandra-BAT spectrum is complex. Our best fit model is the sum of an absorbed 
power-law model and of a reflection component (both having the same photon index of 
2.11^024). The absorbing column density is Niy=3.6l}4 x 10^^ atoms cm"^. The reflection 
component seems to be large R> 1 which is at odds with the absence of the iron Kq, line. 
This source certainly deserves further investigations. 



NGC 3081 is mis-catalogged in SIMBAD as Syl galaxy. In fact the available 6dF spec- 
trum shows clearly that this object is a Sy2 object. We have analized BeppoSax-MECS and 
ASCA data for this source. The best fit is a sum of a black body component, peaking at 
0.58lo'Jf keV, an absorbed power law with column density of 60^3]^ x 10^^ atoms cm"^ and 
photon index 1.91q;q4 ^^"^ line of equivalent width of 24ll}3^eV. 



Table 1. Spectral parameters 



NAME 


RA 
(J2000) 


DEC 
(J2000) 


Type 


r/E[kT]=' 


(lO'^'^ atoms cm^'^) 


MODEL 


INSTR.*" 


oO 1U5.U 


61.9178 


3.6517 


Sy2 


l-66_o.i3 


29.4I4 8 


wabs'^pow 


B, X 


lAXu J0425oo-5711 


66.6021 


-57.1775 


Syl 


-1 +0.028 
-'-•■^^-0.027 





pow 


B, A 


3C 120 


68.2982 


5.3374 


Syl 


1 cn+0.04 /r, 0^+0.026 
1-8U_Q 04/0.2 (•_(, 





wabs*pow+bb 


B, A 


MCG -01-13-025 


72.9205 


-3.8240 


Syl. 2 


1 K+0.48 


< 0.02 (Ij 


pow 


B 


bWil 1 J0505. 7-2348 


"oAbii 


-23.8666 


Sy2 


1 77+0.08 

-■-■"-0.07 


A 0+0.9 


wabs'^pow 


B, X 


CSV 6150 


77.7224 


16.5265 


Syl. 5 


1 nA+0.25 

29.7tl:l 




pow 


B 


4U 0513-40 


78.5146 


-40.0558 


LXB 




brem 


B 


QSO B0513-002 


79.0096 


-0.1332 


Syl 


1 SQ+0.02 r)7+0.02 
l.ec>_o 016/ U-^' -0.02 


< 0.01 


wabs*pow+bb 


B, A 


bWll' 1 J0517.1H-1d33 


79.2839 


16.5605 




r, n+0.23 
^•"-0.26 




pow 


B 


EbU 3d2- (jr Olo 


79.8844 


-32.6720 


Syl. 5 


1 r: + 0.03 
-'-•'J-0.02 


<0.01 


wabs pow 


B,X 


Pictor A 


79.9460 


-45.7557 


Syl 


-1 O+0.015 
-'-■°-0.()14 


r, -, 9 + 0.007 

^•^^-0.02 


wabs pow 


B, A 


ESO 362-G021 


80.6581 


-36.4233 


BL Lac 


-I 7+0.037 
-0.036 


r, 1+0.0197 
^•^-0.0187 


wabs*pow 


B, A, X 


V* TV Col 


82.3541 


-32.7965 


CV-DQ* 


rt+4.6 

24-913.8 




bremss 


B 


V* TW Pic 


83.6470 


-58.0200 


CV 


1 1 c+io.e 




bremss 


B 


LMC X-3 


84.7717 


-64.1148 


HXB 


n r.+0.4 

2-0_o.3 




pow 


B 


LMC A-1 


84.8917 


-69.7210 


HXB 


r, O + 0.22 

^•■^-0.20 




pow 


B 


PSR B0540-69.3 


84.9878 


-69.3230 


Pulsar 


-, r,t:+0.28 
1.85+o2g 




pow 


B 


PKS 0537-286 


84.9953 


-28.7029 


BLAZAR 


l.oo_o 08 


< 0.01 


wabs*pow 


B, A 


r^iYo Uo4o-ozz 


87.7165 




tSLi Lac 


1 Q+0.032 
^•°-0.031 


n 09+0.006 


wabs pow 


rs, A 


NGC 2110 


88.0411 


-7.4554 


Sy2 


^^+0.01 /n .1^+0.02 

l-62Io.oi/0-47Z:o.„2 


A n+0.13 

4-0Io.o7 


wabs*(pow+ga) + bb 


B, A, X 


LEDA 75476 


89.5237 


-38.3799 


Syl 


1 v,i+0.017 / n r,i-+0.08 


9 c + 0.11 
^■''-0.17 


wabs*(pow+ga)+bb 


B, A 


ESO 490- G 26 


100.0031 


-25.8931 


Syl.2 


-1 n/-i+0.05 

1.90^004 


r\ 97+0. 005 


wabs*pow 


B, X 




111 SQf^l 

11 i.oyoi 


9/1 1 (\'iGk 




, cq+0.55 




pow 


D 


V* 441 Pup 


112.1626 


-26.0696 


CV 


12.4l»o« 




bremss 


B 


V* BG CMi 


112.8752 


9.9214 


CV 


qi q+41.2 

ol.o_j^4 2 




bremss 


B 


SWIFT J0732. 5-1331 


113.1328 


-13.5037 


CV 


qq 9+50.1 

■"■^-14.2 




bremss 


B 


SWIFT ,10739.6-3144 


114.9127 


-31.7496 


Sy2d 


1 77+0.51 
-0.43 




pow 


B 


SWIFT J0743.0-2543 


115.7501 


-25.7314 




1 70+0.69 
^■'**-0.56 




pow 


B 


ICR J07597-3842 


119.9822 


-38.7422 


Syl.2 


1 O+0.08 
1-O_0.07 


<0.01 


wabs*pow 


B, X 


UGC 4203 


121.0552 


5.1203 


Sy2 


1 fi8+0 09 / n qi+0.08 
l-68_o.io/ 0.31_o.o5 


i2.5i^:?f 


wabs(pexrav+ga)-|-bb 


B, A, X 


SWIFT J0811.5-I-0937 


122.8750 


9.6214 


XBONG<* 


9 9+2.1 
^■^-0.9 


0'^ 


pow 


B 


SWIFT J0823.4^0457 


125.8271 


-4.9401 


Sy2 


^•°^-0.22 




wabs*pow 


B, X 


Vela PSR 


128.8308 


-45.1771 


PSR 


, 00+0.20 
i.88_o 26 




pow 


B 


FRL 1146 


129.6151 


-35.9976 


Syl 


^■°°-0.31 




pow 


B 



Table 1 — Continued 



NAME 


RA 
(J2000) 


DEC 
(J2000) 


Type 


r/E[kT]=' 


(10^^ atoms cm-^) 


MODEL 


INSTR.'' 




1 on oc^c;ft 


1 O 0/1 fl'/ 

-1Z.z4d / 




1 QC+0.43 


U 


wabs*pow 


J3 


SWFIT J0844.9-3531 


131.2411 


-35.5313 




1 QI+0.46 
^•'^^-0.68 




pow 


B 


SWIFT J0854.7+1502 


133.6828 


15.0371 


Sy2d 




> 0.5" 


pow 


B 


SWIFT J0917.2-6221 


139.112 


-62.359 


Syl 


1 07+0.07 / p, -I .+0.02 


1 qO+0.18 

^■''■^-O.IO 


bb+-wabs*po-w 


B 


Mrk 0704 


139.6505 


16.2987 


Syl.5 


1 36+0-1 
^^•■^°-0.07 


15.0t«-J 


pcfabs* (pow-(-ga) 


B, A, X 


4U 0919-54 


140.0753 


-55.2135 


LXB 


45.1ltf8" 




bremss 


B 


MCG -01-24-012 


140.2134 


-8.0872 


Sy2 


-, 7+0.08 
^•'-0.07 


6 5+0-8 

D.0_o 7 


wabs*pow 


B, X 


NGC 2992 


146.4060 


-14.3007 


Syl.9 


1 24+0-06 
J^-^^-0.05 


17-0-03 


■wabs* (pow-|-ga) 


B, A, X 


ESO 434- G 040 


146.9151 


-30.9388 


Sy2 


1 77+0.006 / n 1 c+0.011 
'-■ ' '-0.07 / ^•^^''-0.016 


, c+0.026 
^•^-0.09 


wabs*(pow+ga)+bb 


B, A, X 


3C 227 


146.9447 


7.4191 


Syl 


2 11+0-14 


3.6 


pexrav+-wa*pow 


B, C 


NGC 3081 


149.8805 


-22.8561 


Sy2 


-I Q+0.02 / ^ c-o+O.lS 
1-9-0.04/ U-08_Q 




wabs*{po-w+ga) +bb 


B, A, S 



^Photon index and/or plasma temperature for the model, specified in column "Model", to fit the data. 

. to 

"^Instrumenst used for spectral analysis are: B = BAT, X = Swift/XRT, A = ASCA, C = Chandra, and S = BeppoSAX. ^ 

Lower limit on absorption estimated through the non detection by ROSAT 
"1 Proposed identification in lRau et al.l ||2007| ). 

"^Order of magnitude of the absorption estimated imposing that the extrapolated source flux match the ROSAT-PSPC count rates. 

^UGC 4203 exhibits transition between reflection-dominated and a trasmission-dominated spectrum. The absorption is estimated in the 
latter case using XRT data (see text for details). 

References. — References: (1) Gallo et al. (2005). 
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Table 2. Spectral parameters for Syl, Sy2, intermediate and all Seyfert AGN. Errors are 

90% confidence level. 



CLASS 


Photon index 


/NDF 


Seyfert 1 


2.23±0.11 


5.4/4 


Seyfert 2 


1.86±0.10 


1.2/4 


Seyfert 1.2-1.5 


1.95±0.11 


4.9/4 


Seyfert All 


2.00±0.07 


2.1/4 



Table 3. Extragalactic sample. 



MA MR 

1 N T\. iV J. Ill 


Tvnp 






Fx 

(lO^"'"^ erg cm~^ '^~^) 


Lx 

flO*''' ere s^-"^") 


T /T a 

^2-lOkeV/^OIII 


(eY) 


IN J/ 

(10^2 atoms cm~2^ 


ref 


3C 105.0 


Sy2 


0.089 


28421 




44.511^ 


581.5±69.2 


nr 


29.4 


1 


lAXG J042556-5711 


Syl 


0.104 


0.78'^ 




55.0^10 




nr 





1 


3C 120 


Syl 


0.0330 


3762 


10.1+?| 


25.4+1-2 


148.8±40.6 


52.3 





1 


MCG -01-13-025 


Syl. 2 


0.015894 


1.15'= 


2.52_ J g 


1 c+o.e 






<0.02 


2 


SWIFT J0505.7-2348 


Sy2 


0.0350 


7.13'= 






345.2±82.2 


nr 


6.3 


1 


QSO B0513-002 


Syl 


0.0327 


0.254 


4.92t°:^ 


I2.3l?;i 




90.8 


0.02 


3 


ESQ 362- G 018 


Syl. 5 


0.0126 


0.58= 




1 7+0.3 
^- ' -0.3 


36.3±3.7 


nr 


<0.01 


1 


Pictor A 


Syl 


0.035 


14045 




5 1+1-5 


113.4±14.5 


nr 


0.12 


1 


ESQ 362-G021 


BL Lac 


0.05534 


2409 


9 7+0.3 
^•'—0 3 


18.815-° 


1284.1±270.3 


nr 


0.1 


1 


PKS 0537-286 


BLAZAR 


3.1 


22000 


2.43l°g 


1.2t°-^ X 105 




nr 


<0.01 


1 


PKS 0548-322 


BL Lac 


0.0690 


383.33 


n 1+0.6 
•^•^-1.0 


37.0l?i°o 




nr 


0.0257 


1 


NGC 2110 


Sy2 


0.007789 


26.92 


27.0l°-g 




1115.8±101.4 


118 


4.0 


1 


LEDA 75476 


Syl 


0.0338 


2.87'= 


q 9+0.6 


8.7111 


393.7±59.6 


144 


2.5 


1 


UGC 4203 


Sy2 


0.01349 


7.67 


4.28+?-^ 


-1 7+0.3 
J^-'-0.6 


214.2±71.4 


747 


12.5 


ito 


SWIFT J0811. 5+0937 


XBONG 


0.282 


268^= 


J^.OO_1.54 


00^+170 
3o4_200 






~0d 


. . 


SWIFT J0823.4-0457 


Sy2 


0.023 


o.ei'^ 


2.78+11 


n ., + 1-1 
■^•■^-1.3 


179.1±71.6 




16.2 


1 


30 206 


QSO 


0.1976 


1194 


2.62+°l 


3001^5^ 









1 


SWIFT J0854.7+1502 


Sy2 


0.0696 






19.7113 






> 0.5'! 


1 


Mrk 0704 


Syl 


0.0292 


0.82^= 


2 21 + ^1 


4 3+2-1 


114.9±27.1 


160 


14.6 


1 


MCG -01-24-012 


Sy2 


0.01964 


2.86'= 


^■"-0.9 


n 7+0.4 
■^•'-0.6 


2914.4±1092.9 


nr 


6.8 


1 


NGC 2992 


Syl.9 


0.00771 


2.03 


3.6+1? 


47+0-13 


0.8±0.1 


520 


0.17 


1 


ESQ 434- G 040 


Sy2 


0.00848 


0.6 


iq 1+0.6 


9 7+0.1 
^•'-0.1 


912.9±21.1 


85.5 


1.5 


1 


3C 227 


Syl 


0.0858 


5462 


2.23l5g 


40.0+1° 


19.3±7.10 


nr 


3.6 


1 


NGC 3081 


Sy2 


0.00798 


0.1'= 


fi S+0-9 
O-»-0.8 


o.96l°:L 


31. Oil. 7 


241 


60 


1 



''The OIII luminosities have been derived in iRau et al.l l|2007h . 

I'Iron line equivalent width. A value of "nr" means that the iron line is statistically not required by the fit. 
'=The radio flux at other wavelengths has been extrapolated to 6 cm assuming fu oc i/~0-5 
Limit on the absorption obtained estrapolating the BAT spectrum to the ROSAT band. 



References. — 



References for the absorption values: (1) this work; (2) Gallo ct al. (2005); (3) Lutz et al. (2004). 



-30- 



Table 4. Comparison with previous results 



Instrument 


Ref. 


Energy 


AGN density* 


BAT density (this work)^ 






keV 


10-2 deg-2 


10-2 deg-2 


INTEGRAL-ISGRI 


1 


20 - 


40 


0.48±0.08 


0.41ib0.08 


INTEGRAL-ISGRI 


2 


100- 


150 


0.18±0.006 


0.17±0.034 


HEAO-1 A2 


3 


2 - 


10 


1.2±0.2 


1.6 ±0.32 


RXTE PGA 


4 


8 - 


20 


0.56±0.06 


0.65±0.13 


XMM 


5 


0.5 


- 2 


O.liO.Ol 


1.3 ±0.26 


XMM 


5 


2 - 


10 


0.95±0.06 


1.6 ±0.32 


XMM 


5 


5 - 


10 


0.63±0.4 


0.4±0.08 



*AGN densities from different surveys above 10 erg cm ^ s ^ (in the respective 
bands) . 

^The BAT AGN density was converted to the native energy band of the measure- 
ment we are comparing it with. 



References. — (l)Beckmann et al. 2006; (2) Bazzano et al. 2006; (3) Piccinotti et 
al. 1982; (4) Revnivtsev et al. 2004; (5) Gappelluti et al. 2007 
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Fig. 1. — Plot of HRl and HR2 hardness ratios. The sohd hne is the locus for sources with 
unabsorbed power law spectra with photon indices from 1.0 to 3.0 while the long dashed 
line shows the location of Compton-thick AGN with the same range of photon indices. The 
dashed line shows the location of objects with a thermal bremsstrahlung spectrum with 
temperatures in the range 5-50 keV. In the upper left corner the typical ±lcr error for a 5 a 
source is shown. 
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Fig. 2. — Stacked spectrum of all AGN reported in Tab. [3] excluding the Blazars. The 
dashed line is the best power-law fit to the data (photon index 2.0±0.07.) 
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Fig. 3. — V/Vmax as a function of detection threshold for the sample of extragalactic 
sources. The dashed line is the expected value (0.5) for a complete sample in an homogeneous 
distribution. The solid line shows the mean test value for S/N> 4.5o". 
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Fig. 4. — Extragalactic cumulative source count distribution in the 14-170 keV band. The 
sohd hne is the best fit described in the text. 
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Fig. 5. — Luminosity, in the 14-170 keV band, vs. intrinsic column density for the extra- 
galactic sample. The blazars are highhghted with a triangle. 
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Fig. 6. — Crab rates in the 14-22 keV band as a function of the tangent of the off-axis angle. 
When the Crab is 50° off-axis the detected count rate is ~30% lower than the on-axis count 
rate. The solid line is a polynomial fit to the rates. 
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Fig. 7. — Folded spectra and best fit models as described in the text. From left to right and 
up to bottom the spectra are for: 3C105.0, lAXG J042556-5711, 3C120, MCG-01-13-025, 
SWIFT J0505. 7-2348 and CSV6150. 
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Fig. 8. — (for online version) Folded spectra and best fit models as described in the 
text. From left to right and up to bottom the spectra are for: 4U 0513-40, QSO B0513-002, 
SWIFT J0517.1+1633, ESO 362- G 018, Pictor A and ESO 362- G 021. 
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Fig. 9. — (for online version) Folded spectra and best fit models as described in the text. 
From left to right and up to bottom the spectra are for: TV Col, TW Pic, LMC X-3, LMC 
X-1, PSR B0540-69.3 and PKS 0537-286. 
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Fig. 10. — (for online version) Folded spectra and best fit models as described in the text. 
From left to right and up to bottom the spectra are for: PKS 0548-322, NGC 2110, LEDA 
75476, ESO 490- G 26, J0727.5-2406 and V441 Pup 
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Fig. 11. — (for online version) Folded spectra and best fit models as described in the 
text. From left to right and up to bottom the spectra are for: BG CMi, J0732. 5-1331, 
J0739.6-3144, J0743.0-2543, IGR K07597-3842 and UGC 4203. 
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Fig. 12. — (for online version) Folded spectra and best fit models as described in the text. 
From left to right and up to bottom the spectra are for: J0811. 5+0937, J0823.4-0457, VELA 
PSR, FRL 1146, 3C 206 and J0844.9-3531. 



-41 - 



SWIFT J0S54.7+ 1502 



SWIFT J0917.2-6221 



channel energy (keV) 



channel energy (keV) 



MCG -01-24-012 



I i 
1 ^ 

I h 

I 6 



channel energy (keV) 





> 6 

I I 

I h 

1 i 

i 6 

o 




channel energy (keV) 



4U 0919-54 



g 6 



channel energy (keV) 



I 6 

a 6 



channel energy (keV) 







- 1 





Fig. 13. — (for online version) Folded spectra and best fit models as described in the text. 
From left to right and up to bottom the spectra are for: J0854. 7+1502, J0917.2-6221, Mrk 
704, 4U 0919-54, MCG -01-24-012 and NGC 2992. 
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Fig. 14. — (for online version) Folded spectra and best fit models as described in the text. 
From left to right and up to bottom the spectra are for: ESO 434 -G 040, 3C 227 and NGC 
3081. 



